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ABSTRACT 

A new method to constrain the local non-linear coupling parameter /jvl based on a 
fast wavelet decomposition is presented. Using a multircsolution wavelet adapted to 
the HEALPix pixelization, we have developed a method that is ~ 10 2 times faster than 
previous estimators based on isotropic wavelets and ~ 10 3 faster than the KSW bis- 
pectrum estimator, at the resolution of the Wilkinson Microwave Anisotropy Probe 
(WMAP) data. The method has been applied to the WMAP 7-yr V+W combined 
map, imposing constraints on /jvx of —69 < $nl < 65 at the 95 per cent CL. This 
result has been obtained after correcting for the contribution of the residual point 
sources which has been estimated to be A/tvl = 7 ± 6. In addition, a Gaussianity 
analysis of the data has been carried out using the third order moments of the wavelet 
coefficients, finding consistency with Gaussianity. Although the constrainsts imposed 
on /tvl are less stringent than those found with optimal estimators, we believe that 
a very fast method, as the one proposed in this work, can be very useful, especially 
bearing in mind the large amount of data that will be provided by future experiments, 
such as the Planck satellite. Moreover, the localisation of wavelets allows one to carry 
out analyses on different regions of the sky. As an application, we have se parately anal- 
ysed t he two hemispheres defined by the dipolar modulation proposed bv lHoftuft et al.l 
(12009ft . We do not find any significant asymmetry regarding the estimated value of 
Jnl in those hemispheres. 
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1 INTRODUCTION 

The cosmic microwave background (CMB) is one of the pil- 
lars that provide support to the Big Bang theory. The fluc- 
tuations of the CMB naturally arise in an inflationary sce- 
nario. The understanding of this very early stage of the his- 
tory of the Universe is a challenging issue for the scientific 
community due to the implications on large scale structure 
formation and fundamental particle physics at high ener- 
gies. A large number of inflationary models have bee n pro- 
posed in the literature (for an overview see for instance lLvthl 
2008) but the task of testing such scenarios is not trivial, and 
there is the need of new experiments and powerful statistical 
tools to discriminate among them. In this sense, the statis- 
tical properties of the CMB temperature anisotropics are a 
source of information about the processes that have gener- 
ated the primordial fluctuations. In particular, the standard, 
slow roll, single field inflationary model predicts a nearly 
Gaussian distribution of the CMB temperature anisotropies, 
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while alternative models may introduce a certain level of 
non-Gaussianity in the CMB. A convenient parametriza- 
tion valid for a large set of non-standard inflationary mod- 
els which includes the quadratic corrections of the pri- 
mordial curvature perturbation is dSalopek fc Bondlli 990: 
Gangui et ail 1 1994 IVerde et al.1 l200d : iKomatsu fc Spergel 
200lf ): 



(r) = Mr) + fNL [<pl-<<pl >] 



(1) 



where <j>L are Gaussian linear perturbations and }nl char- 
acterises the amplitude of the non-linear contribution in 
real space. This local for m appears in non-standard multi- 
field inflationary models l|Babich et al.|[20o3 : IKomatsu et al.l 
2009). For a complet e review on non-Gau ssianity due to in- 
flationary models see iBartolo et all (|2004l )- 

In addition to in- 
flationary models, there are other alternative scenarios that 
can be constrained, such as the ekpyrotic mo del where a 
negative value of Jnl is expected (|Lehnerdl2oioh . Moreover, 
there are other processes that can introduce deviations from 
Gaussianity in the third order moments (as foreground con- 
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tamination, non-linear gravitational effects, topological de- 
fects, etc). 

Since the quadratic parametrization was proposed, 
an important effort has been made to set observational 
constraints on local /jvl with a wide variety of meth- 
ods including the bispectrum (jYadav fc Wandeltl 120081 ; 
ISmith et all 12009"; IKomatsu et al.l |2010|), wavelet - based 



methods JCavon et al.1 I200&1; iMukheriee fc Wand |2004 



Curto et al.1 l2009al lbl; iPietrobon et al.l 120091; IRudiord et al] 
20091 ), Minkowski functionals dHikage et al.1 120081 ) or the 
N-pdf (|Vielva fc Sanj|2009l . l2010l ). Most of these works find 
that the data are compatible with /jvl =0, but the con- 
straints are not yet sufficiently tight to discriminate among 
a large set of inflationary model s. The best current limit 
is given by IKomatsu et al.1 (|2010t ) and is -10 < /jvl < 74 
at the 95 per cent confidence level. These constraints have 
been obtained with a bispectrum estim ator, which i s com - 
putationally very demanding. However, ICurto et al.l (|2010l ) 
have recently shown that an estimator based on the SMHW 
can provide constraints on /jvl as stringent as the optimal 
estimator, the bispectrum, while reducing considerably 
the CPU time. With the arrival of new data from high 
resolution experimen ts such as the ESA Planck satellit^] 
l|Tauber et al.l |2010|) . it becomes even more important 
the availability of even faster and simpler methods. With 
this aim we present the application for CMB of a wavelet 
adapted to the HEALPix pixeliz ation similar to the tool 
proposed bv lShahram et al. (2007). 

The paper is organised as follows. In Section [2] we in- 
troduce the HEALPix wavelet decomposition. In Section [3] 
the method to constrain the /jvl parameter as well as the 
proposed Gaussianity test are described. In Section [5] we 
present the results of the application of this technique to the 
WMAP 7-yr data. Finally, our conclusions are summarised 
in Section [S] 



2 THE HEALPIX WAVELET 

A large set of different wavelets have been used in 
the astrophysics literature. In particular, different spher- 
ical wavelets have been applied to CMB Gaussianity 



last decade, 



analysis during the 

ical Ha ar wavelet (SHW, iBarreiro et al, 
SMHW |Martmez-Gonzalez et alj |2002|; IVielva et al.|[200i; 



including the spher- 
|2000|) the 



Mukheriee fc Wand I200I ; ICruz et al l 120051; ICurto et al.1 
2009bl ). elliptical SMW H iMcEwen et al.1 I2OO5I) direc - 
tional spherica l wavelets (IMcEwen et al.ll2005l. [2 006, 2008) 
and needlets dPietrobon et al.l 120091 ; IRudiord et al.1 12009 ; 
ICabella et a"i1l2010l). For a review on wavelet applications 
to cosmology see lMcEwen et ai1l2007l . 

In this work, we will present an application using the 
so-called HEALPix w avelet (HW). In the previous work of 
l|Shahram et al.1 120071 ) a linear operator is applied to the 
HW to obtain wavelet coefficients corresponding to vertical, 
horizontal and diagonal orientations. This operation leads 
to a wavelet coefficients without redundancy, obtaining a 
number of wavelet coefficients (details plus approximation) 
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equal to the number of original pixels. However, we have 
kept the HW with its intrinsic redundancy for three main 
reasons: first, to improve the computational time, second 
to obtain a wavelet decomposition as isotropic as the HW 
allows 0, and third, because, as it is shown later, redun- 
dancy helps to improve the sensitivity in the detection of 
/jvl- Similarly to the SHW, the HW is a discrete, orthog- 
onal wavelet, adapted to a hiera rchical pixelization (such 
as HEALPbPI. lG6rski et alj|2005l ). whereas the SMHW is 
a continuous, non-orthogonal wavelet and does not have a 
hierarchical decomposition structure. The HW presents an 
optimal space localization, while the scale localization is not 
as good as that of the SMHW. It is important to point out 
that, in the case of the SMHW and of needlets, a transfor- 
mation of the data into spherical harmonic space is required. 
However, this is not the case for the HW and, therefore, the 
computational cost is significantly reduced. 
The resolution of a HEALPix map is characterised by the 
TVsidc parameter, such that the number of pixels in which 
the sphere is divided corresponds to TV = 12/V s 2 ide . iV 3 ide 
can only take powers of two as values. The HW decomposes 
the temperature map at resolution J, where 7V s id c = 2 , in 
wavelet coefficient maps at all the allowed HEALPix resolu- 
tions down to the lower considered resolution jo - The wavelet 
functions are defined as follows: 



®l.j,k( x ) = (fj+i,k 1 (x) - 

^2,j,k(x) = (p j + 1M (x) - 

where ip(x)j t k is the scaling function 



<PjA x ) 
4 

fj,k{x) 
4 

<Pj,k( x ) 
4 

<Pj,k( X ) 



(2) 



tp(x) 



j,k 



1 if X £ P jtk 
otherwise , 



(3) 



and Pj t k is the pixel at position k at resolution j, which 
at the next higher resolution is divided into four daughter 
pixels Pj +liko , Pj + i jftl , P j+1 , k2 , P j+ i tk3 . 

The wavelet decomposition of a temperature map can 
be written in terms of the basis functions and a set of 
coefficients: 



AT, 



k=0 
J 3 

+ ^ ] ^ ] ^ ] 'Jm,j,k^m,j,k{Xi) , 
j~jo m=0 fe=0 



(4) 



where Nj is the number of pixels at resolution j. \j tk and 
7m,j,fe are the approximation and detail coefficients respec- 
tively. From a practical point of view, to perform the decom- 
position, we start with the original resolution, i.e. j — J. 



2 It is worth mentioning that HW detail coefficients help to 
highlight the isotropy properties of the field as compared to the 
directional oriented details of the SHW. This is important because 
the local non-Gaussianities are expected to be isotropic. 

3 http://healpix.jpl.nasa.gov/ 
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At this resolution, the approximation coefficients Aj^ cor- 
respond to the pixels of the original tempature map. The 
approximation coefficients at the next resolution are simply 
obtained by degrading the map to the inmediatly lower res- 
olution (i.e., by averaging the corresponding four daughter 
pixels): 



A 



! 3 



(•>) 



spectrum that best fits the WMAP-7yr data, and that was 
used for the Gaussian simulations (the original simulations 
were obtained using the WMAP 5-yr power spectrum). 
Again, we construct the maps for the V and W receivers, ap- 
plying the corresponding beam and pixel transfer functions 
and adding the appropriate level of noise. Finally a single 
V+W combined map is constructed for each non-Gaussian 
simulation. 



On the other hand, the detail coefficients at resolution j + 1 
are simply obtained by subtracting the approximation at res- 
olution j from the approximation at resolution j + 1. Thus, 
the detail coefficients are defined as: 



7o,j,fc 


= A 3 - 


n,k 


~ 4Aj,* 


(6) 


Ji,j,k 


= V 


hi,*!! 






72,j,fc 


= V 


rl,k 2 






73,j,fc 


= V 


n,k 3 


- 4A,, fc 





A schematic diagram of how to obtain the approximation 
and detail coefficients is given in Fig. [T] 



3 METHODOLOGY 

The main purpose of this work is to constrain the parame- 
ter f NL defined in Eq. (JTJ) using the WMAP-7yr datfJE For 
this analysis, we only consider the (foreground reduced) V 
and W channels, since they are less afected by foreground 
contamination. A single CMB map is obtained through a 
noise-weighted linear combination of the V and W receivers. 
The KQ75 mask (which covers around a 29 per cent of the 
sky) is applied subsequently. 

In order to calibrate our estimator, we need both Gaus- 
sian and non-Gaussian simulations. To generate the Gaus- 
sian simulations, we compute the power spectrum that best 
fits the W MAP-7yr data a c cordin gly to the parameters es- 
timated bv lKomatsu et all (|2010|). For this purpose we use 
the On-line tool CAMB |Lewis et al.ll200(J ). We then apply 
the corresponding beam and pixel functions to simulate the 
data at each of the considered receivers (2 for V and 4 for 
W). A Gaussian noise realisation is subsequently added to 
the CMB maps with a variance per pixel given by ,J ° , 
where oq is the detector sensitivity of each of the receivers 
and Nobs is the number of observations at each pixel. Finally 
the six maps are combined in the same way as the data. 

Regarding the non-Gaussian simulations, we have used 
the 1000 simulations generated bv lElsner fc Wandeltl (|2009l 'l 
that are publicly available^. The previous authors provide 
the harmonic coefficients for the Gaussian and non-Gaussian 
parts of the simulation. A non-Gaussian simulation with a 
given value of /jvl is then constructed as: 



Oil, 



AG) 



, t (JVG) 



(7) 



where we have normalized and a;^ G '0 to the power 



4 The data are available at the LAMBDA web page: 
http: / /lambda. gsfc.nasa.gov/ 

b http:/ /planck. mpa-garching.mpg.de/cmb/fnl-simulations/ 

6 The amplitude of the a,^* has been corrected by a factor of | 

as indicated by the authors 



3.1 Cubic statistics 

In this section we define the third order moments of the 
wavelet coefficients that are used to constrain Jnl- Sim- 
ilar statistics have been used in oth er previous works 
|Curto et al.ll2009bl : iRudiord et al.ll2009l ). 
We perform the wavelet decomposition of the considered 
map starting at resolution A^idc = 512 ( J = 9) and down to 
iV s idc = 2 (jo = 1) (when using a higher value of jo we are 
losing efficiency whereas for jo = the results are not sig- 
nificantly improved while the computational time increases 
by a 30%). We obtain 8 detail maps and 1 approximation 
map. In addition, we also include in the analysis the origi- 
nal map and the 8 intermediate approximation maps (which 
are obtained during the wavelet decomposition to construct 
the detail coefficients). Although, in principle, these addi- 
tional approximation maps contain redundant information, 
they seem to provide additional information regarding the 
third order statistics, since a larger number of third order 
combinations can be constructed. In fact, we have tested 
that with the inclusion of the approximation maps and the 
original map, the results are improved by a 30%. Therefore, 
we have a total of 18 maps for each analysed signal. The 
statistics are constructed as the third order moments of all 
the possible combinations of these 18 maps, where the coef- 
ficients are weighted to take into account the presence of a 
mask. In order to calculate these weights, one performs the 
wavelet decomposition of the considered mask (that has ze- 
ros in the masked pixels and ones in the rest). The wavelet 
coefficients of the mask at each detail and approximation 
scale are used to construct the weight Wj(i) of the coeffi- 
cient at position i at resolution j. This makes sense if one 
bears in mind how the wavelet decomposition is performed. 
For instance, to construct the approximation map at resolu- 
tion J — 1 at a given position i, the four daughter pixels at 
resolution J have to be averaged. If the four pixels are un- 
masked, this corresponds to a weight of 1 in the original map 
and also in the approximation map at position i. However, 
if one of the orignal pixels is masked, this pixel would have 
a weight of zero, and the average would be done only over 
three pixels. Thus the weight of the corresponding approx- 
imation coefficient would be 3/4. Therefore, this weighting 
scheme takes into account the fact that different coefficients 
contain different amount of information, depending on the 
considered mask. Also, contrary to the case of other wavelet 
estimators, this means that the mask does not need to be ex- 
tended but, in fact, it is reduced when increasing the scale. 
This is due to the fact that a larger pixel is kept for the 
analysis, with the appropriate weight, if at least one of the 
daughter pixels was unmasked. 
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Figure 1. Diagram of the construction of the approximation and detail coefficients. Approximation coefficients are computed as the 
average of the four daughter pixels. Detail coefficients are computed as the subtraction of that average from the original pixels and are 
represented by d. 



The third order statistics are then defined as: 



3.2 Gaussianity test and }nl constraints 



JVj-l 



E 



UjtJkCri 



(8) 



where 



= 2/i. 



fij and yi t j are the wavelet coefficients 



maps at position i at resolution j. Note that j goes from jo 
to J, k goes from j to J and I goes from k to J. fij and 
ex, are the weighted mean and the dispersion for the map at 
resolution j. Wjki(i) is the weight associated to the wavelet 
coefficients at position i and scales j, k, I and is given by: 



Wjki(i) = \/wj(i)w k (i)wi(i) 



(9) 



We first perform a Gaussianity test in order to probe 
whether the data is compatible with Gaussianity using the 
X 2 estimator: 



(10) 



»ij=i 



where Vi is the vector of the third order statistics computed 
from the considered map (to simplify notation, hereinafter 
we define Mj = S'm,V2 = S113,...). (vi) and dj are the 
mean and covariance matrix of the statistics obtained with 
10000 Gaussian simulations. To perform the Gaussianity 
test, the value of the x 2 ls computed for the WMAP data, 
and compared to the distribution of the estimator obtained 
from an additional set of 1000 Gaussian simulations. 



Note that some of these statistics are redundant (lin- 
early dependent between them), so we restrict our analysis 
to the set of non-rendudant statistics, which gives a total of 
fistat =232 quantities. 

The process for computing these statistics requires ~ 
N x rz s tat number of operations, where N is the number 
of pixels and n s tat the number of statistics computed. This 

number is significantly lower than that of the full bispectrum 

5 

that nee ds N 2 operations. Usi ng the KSW algorithm pre- 
sented in lKomatsu et all (|2005t ) the number of operations is 

3 

reduced to ~ rN 2 , where r is the number of sampling points 
(of the order of 100). On the other hand, the SMHW scales 
as ~ n s N 2 , where n s is the number of scales involved (of 
the order of 10). Thus, at WMAP resolution (JV sid e = 512 
and ~ 3 x 10 6 ) we have that the method presented in this 
work is 10 2 times faster than the SMHW, 10 3 times faster 
than KSW bispectrum estimator and 10 7 faster than the 
general bispectrum estimator. 



The second analysis that has been performed is the es- 
timation of /nl from the data. As the wavelet decompo- 
sition is linear, we can obtain the wavelet coefficients from 
the Gaussian and non-Gaussian parts separately. Thus, the 
wavelet coefficients for a given value of /nl are given by 

(G) . r (NG) / 1n \ 

Vi =Vi +fNLVl ■ (11) 

Taking into account that y< NG > are around 4 orders of mag- 
nitude smaller than 2/ , when we compute < y 3 > the NG 
high-order terms can be neglected and it can be shown that 
/jvl is proportional to the wavelet estimators, as it is also 
the case for other statistics (such as the bispectrum): 

Vi = CLifNL , (12) 

where a; can be computed from simulations with a simple 
linear regression. 

In order to estimate the }nl parameter, we perform a 
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X 2 minimisation. In particular, x 2 (/jvl) is defined as follows 

Tlstat 

= ( Vi ~ (""^fNL^r/ifNL)^ - (vj) fNL ) , (13) 

where (vi)f NL is the mean of the statistics for a given 
value of /jvl obtained from the 1000 non-Gaussian simula- 
tions and Cy(/jvx) is the corresponding covariance matrix. 
For /jvl << 1000 is reasonable to use the approximation 
Cij(/jvz,) — Cij, where dj is the covariance matrix for the 
Gaussian case. 

Error bars on the parameter estimation at different con- 
fidence levels are found using the Gaussian simulations. We 
also compute the minimum variance in a semi-analytical 
manner. It is well known that the diagonal of the inverse 
of the Fisher matrix provides an estimation of the variance 
of the parameters. In order to estimate the Fisher matrix, we 
approximate the distribution of the statistics by a Gaussian. 
Using this approximation and taking into account Eqs. (|12p 
and (|13[) , the variance from the Fisher Matrix can be written 
as: 



E 



diC i: j dj 



In practice, the distribution of the statistics do not follow a 
perfect Gaussian distribution. Therefore, this variance can 
be seen as a lower limit to the true underlying variance. 



4 RESULTS 

In this section, we present the analysis of the WMAP-7yr 
V+W combined map. On the one hand, we analyse the 
compatibility of the data with Gaussianity using the cu- 
bic statistics defined in Eq. (JH| and the estimator presented 
in Eq. (|10|l . On the other hand, we compute the best-fit 
/jvl parameter from the data by minimizing Eq. (|13[) . Error 
bars are set using simulations. In addition, we also present 
a study of the contribution of the point sources to the es- 
timated /ml value and of the variation of / jvl estimated 
from two independent hemispheres (defined bv lHoftuft et al.l 
120091) . 



4.1 Gaussianity test 

As explained in Section 13.11 we have considered a total of 
232 cubic statistics, constructed from 18 maps at 9 differ- 
ent scales. Fig. [2] shows the value of Vi for the WMAP 7-yr 
V+W data, after applying the KQ75 mask. The plot does 
not show any obvious deviation from Gaussianity. To fur- 
ther study the consistency of the data with Gaussianity, we 
also perform the \ 2 t es t defined in Eq. (|10[) . From 1000 
Gaussian simulations, we estimate the distribution of this 
quantity, finding a mean value of (x 2 ) = 233, very close to 
the number of degrees of freedom (232). The value of the 
dispersion is 69, larger than expected for a \ 2 distribution 
with the considered degrees of freedom. However, this may 
be explained by the fact that the distribution of the dif- 
ferent statistics are not purely Gaussian. For the WMAP 
data, we find Xdata = 434 with a cumulative probability of 




Figure 2. The cubic statistics Vi from WMAP-7y V+W data are 
shown. Shadow areas correspond to the 68, 95 and 99 per cent 
confidence levels of the distribution obtained from 1000 Gaussian 
simulations. The statistics have been plotted from lower to higher 
variance. 




P(x 2 



Xdata 



0.96. Although the result indicates that 



Figure 3. Mean values of the cubic statistics m from 1000 non- 
Gaussian simulations with /jvl = 0, ±100, ±300. 



the WMAP data is some how in the tail of the distribution, 
the x 2 value is not large enough to claim a deviation from 
Gaussianity. 

4.2 Constraints on the /jvl parameter 

We have also performed an estimation of the non-linear 
parameter /jvl- As already mentioned, for this analy- 
sis we h ave used the 1000 non-G aussian simulations pro- 
vided by El sncr fc Wandeltl l|2009h . In Fig. El the mean of 
the cubic statistics derived from simulations with /jvl = 
0, ±100, ±300 is presented. It can be seen that, as stated in 
Eq. (|12[1 . the statistics are proportional to the value of /jvl. 

After obtaining the cubic statistics for the WMAP-7yr 
data and minimising the x 2 given by Eq. (|13|l . we estimate 
that the best-fit value of /jvl is 6. Using Gaussian simu- 
lations, we find that the contraints for the parameter are 
—28 < /jvl < 40 at the 68 per cent confidence level and 
—62 < /jvl < 72 at the 95 per cent confidence level. It is 
also interesting to point out the agreement between the dis- 
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persion computed semi-analytically through the Fisher ma- 
trix (Eq. I14|) and that obtained from Gaussian simulations, 
which are both estimated to be around 34. 

Although the constraints provided by the HW are less 
stringent than those found with optimal estimators (such as 
the bispectrum or the SMHW), they are similar or even bet- 
ter than those obtained by other mehtods such as needlets 
l|Pietrobon et alj|2009l; IRudiord et al]|2009h . the Minkowski 
funct ionals l)Hikage et al.ll2008t ) or the N-pdf (|Vielva fc Sand 
2009). Moreover, as already pointed out, our estimator is sig- 
nificantly faster than all the previously mentioned methods, 
providing a very valuable tool, especially for future high res- 
olution experiments such as Planck. It is also interesting to 
point out that we find a more symmetric co nstraint around 
zero t ha n those obtained, f or instance, by iKomatsu et al.l 
(|2010h or lCurto et ail (|201Ch . 

In order to study further the robustness of our estima- 
tor, we have performed some additional tests. In particular, 
we have estimated the mean value and dispersion of the 
best-fit value of /jvl from simulations with different values 
of /jvl- The left panel of Fig. 0] shows the histograms of 
the estimated values of /jvl obtained from simulations with 
/jvl =-40,0,20,60. To carry out these tests, we have used 500 
of the 1000 non-Gaussian simulations to estimate the mean 
value of the third order statistics (vi)f NL and the remaining 
500 simulations to obtain estimates of /jvl and construct 
the histograms. The mean values and dispersions of /jvl are 
given in the corresponding panels. In particular, we see that 
the method is unbiased, since the mean value of the esti- 
mated /jvl is very close to the true underlying value for all 
the considered cases. In addition, we also plot in the right 
panel of Fig. [4] how the dispersion of the estimator varies 
as a function of /jvl- The standard procedure to estimate 
this dispersion is to use Gaussian simulations but, as seen 
in the plot, this gives a minimum in the estimated value of 
a (/jvl). However, for small values of /jvl, such as the ones 
found in this paper, the variation is small and therefore one 
can safely use the value of the dispersion estimated for the 
Gaussian case. 

Finally, we have repeated the same analy sis using a set 
of 300 non-Gaussian simulations generated bv lLiguori et al.l 
(2007), finding very similar constraints on /jvl- 



4.3 Point source contribution 

The background of unresolved point sources may introduce 
a bias in the estimation of /jvl- In order to correct this 
bias, we have studied the contribution to /jvl given by a 
point source background that is added to the Gaussian sim- 
ulations. For that purpose, we have p roduced point sourc e 
simulations following the procedure of ICurto et al.l (|2009aT ). 
In particular, point sources maps are simulated accord ing to 
the density distribution given bv lde Zotti et all (|2005T ) in a 
range of intensities between I m in — lmjy and I m ax = Uy. 
These maps are then convolved with the corresponding beam 
and pixel functions and added to the simulations contain- 
ing Gaussian CMB plus noise. The estimated value of /jvl 
when point sources are present is then compared to the one 
obtained from simulations containing only CMB and noise, 
finding a difference of A/jvl = 7±6. Fig. [S] shows the effect 
that point sources have on the V{ statistics. As one would ex- 
pect, they mainly affect the statistics involving small scales, 




Figure 5. Mean values of the Vi statistics obtained from 1000 
Gaussian simulations with and without point sources. Diamonds 
represent CMB plus noise simulations, while crosses correspond 
to simulations including also the point sources. To improve the 
visualization, the statistics have been normalised to unit disper- 
sion. 



that correspond to the ones with a lower value of k in the 
figure. Taking into account this result our final constraint 
on /jvl for the WMAP-7yr data is -69 < /jvl < 65 at the 
95 per cent confidence level. 



4.4 Local study of /jvl 

Finally, we have analysed the data considering two in- 
dependent hemispheres. In particular, we have consid- 
ered the hemispher es associated to th e dipolar modu- 
lation proposed by Irloftuft et al.l (|2009l ) where the pre- 
ferred direction is pointing towards the Galactic coordi- 
nates (l 1 b)=(224° J -22 t> ). We have estimated the best-fit 
value and constraints on /jvl for the WMAP 7-yr data 
in both hemispheres, following the same procedure as for 
the full-sky. After correcting the point source contribu- 
tion, the constraint found for the northern hemisphere is 
—73 < /jvl < 119 while for the southern hemisphere we 
have —137 < /jvl < 62, both at the 95 per cent confidence 
level. Therefore, as it was the case for the full-sky, both hemi- 
spheres are consistent with Gaussianity (i.e, /jvl = 0). We 
have also tested that the results from the hemispheres are 
consistent between them. In particular, we have obtained the 
mean difference and dispersion between the /jvl estimates 
at each hemisphere for Gaussian simulations, finding values 
of (A/jvl) =-4 and a (A/jvl) =71. For the WMAP data, we 
have A/jvl = 67, which is perfectly consistent with the val- 
ues expected from Gaussian simulations. Therefore, we do 
not find any assymmetry for the considered hemispheres. 
These results are in agreement w i th th e an alysis based on 
needle ts made bv lPietrobon et ail (|2010h and lRudiord et ail 
(2010) for the WMAP-5yr data, where several divisions of 
the CMB map are studied w ithout finding a sign ificant 
asymmetry. In a recent work, IVielva fc Sand (|2010i ) have 
found an asymmetry in the same hemispheres studied in 
this paper on the estimation of the /jvl using the N-pdf. 
The disagreeme nt may be caused by th e dif ferences on the 
methods. While iPietrobon et al.l (|2010T > and IRudiord et al.l 
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Figure 4. The left part of the figure shows the histograms of the estimated /jvl from simulations with values of /jvl = —40,0,20,60. 
The mean value and dispersion of /jvl for each considered case is indicated in the corresponding panel. The right panel shows the 
behaviour of 0"(/jvl) when estimated from simulations with different values of /jvl- 



(2010) h ave worked with the s ame resolution as we did (6.9 
arcmirO. lvIelva fc San J (|201fjh focused on scales around 2°. 
Also, the non-Gaussian model used by the former works is 
the same as the one used in this paper, whereas the model 
of the latter stands on the Sachs- Wolfe regime. 



5 CONCLUSIONS 

We have presented a new methodology to analyse the Gaus- 
sianity of the CMB and to constrain the /jvl parameter 
using the so-called HEALPix wavelet. To our knowledge, 
the developed /jvl estimator is the fastest method that has 
been proposed up to date. In particular, for WMAP resolu- 
tion (Aside = 512), it is ~ 10 2 times faster than the SMHW, 
~ 10 3 times faster than the KSW bispectrum and 10 7 times 
faster than the general bispectrum estimator. Moreover, al- 
though the constraints imposed by our method are not as 
stringent as those of the optimal estimators (based on the 
bispectrum or on the SMHW) , they are very similar or even 
better than those proposed by alternative methods, such as 
needlets, Minkowski functionals or the N-pdf. 

The method, which is based on the calculation of the 
third-order moments of the wavelet coefficient maps, has 
been applied to the WMAP-7yr V+W combined map. On 
the one hand, we have performed a x 2 test to study the 
Gaussianity of the CMB, finding consistency with the Gaus- 
sian hypothesis. On the other hand, we have constrained the 
value of the local /jvl parameter to be —69 < /jvl < 65 at 
the 95 per cent confidence level, after correcting for the point 
source contribution. In addition, the HEALPix wavelet gives 
the possibility of performing local studies of Gaussianity in 
the CMB map. In particular, we have analysed two inde- 
pendent he mispheres assoc i ated to the dipolar modulation 
proposed bv lHoftuft et alj (|2009T ). In this study, we do not 
find any significant asymmetry on the /jvl estimates for the 
two hemispheres of the WMAP data. The constraints for the 
northern and southern hemispheres are —73 < /jvl < 119 
and —137 < /jvl < 62, respectively, at the 95 per cent con- 
fidence level. 
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